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1 Introduction 



The scattering of charged particles by an infinitely long straight solenoid that encloses a mag- 
netic flux, known as the Aharonov-Bohm effect, is of paramount interest in quantum physics Q. 
Even though the region containing the magnetic field is inaccessible to the particles, the magnetic 
flux inside the solenoid affects their propagation. The observed interference pattern cannot be ex- 
plained within classical physics; it is a purely quantal effect, without classical correspondence. This 
fact singularizes the Aharonov-Bohm (AB) effect from other important processes, like Rutherford 
scattering. 

The AB effect has been analyzed by many different approaches. We shall restri ct ourselves 
to the simplest situation of a straight filiform solenoid with constant magnetic flux. The object of 
interest is the wave function of the scattered particles and the ensuing scattering amplitude. There 
is some freedom in choosing the initial wave function, but we shall consider onl y two possibilities: 
a plane wave or a well-localized wave packet. In the first case we can drop the time dependence 
of the wave functions without losing any information; this situation will be referred to as "time- 
independent scattering". In the second case we shall speak of "time-dependent scattering". Most 
of the early analysis of the AB scattering can be considered time-independent. 

The scattering amplitude was first calculated by Aharonov and Bohm |jj by means of a de- 
composition of the wave function in positive, negative and zero angular momentum components, 
which are calculated separately. The long-distance form of the scattered wave far from the forward 
direction defines the scattering amplitude that bears these authors' names. They did not determine 
the wave function along the forward direction, where the scattering amplitude see ms to diverge. 

A different procedure is due to Berry § , who obtained an exact single- valued wave function by 
suming the wave function endowed with the Dira c phase factor over all contributions ("whirling 
waves") representing differe nt windings round the flux line. It is also possible to obtain the same 
result b y including the geometrical phase factor (Berry's phase) instead of summing ove r whirling 
waves ||. In these works, which have provided a thorough understanding of the AB effect, the 
wave function is presented as a series of Bessel functions which is not further elaborated. 

A more explicit calculation of the AB wave function was supplied by Takabayashi . An integral 
representation for the Bessel functions enabled this author to evaluate the asymptotic limit of the 
wave function far from the forward direction, thus reproducing the AB scattering amplitude, and 
the leading term of the wave function near the forward direction. Both the scattered and the 
transmitted wave were shown to be discontinues along the forward direction, in such a way that 
the total wave function was continuous. 

In the context of a more general investigation, Jackiw ||] found an integral representation for 
the total AB wave function similar to the one found in 2+1 dimensional gravity [||. His integral 
representation allows for a natural decomposition of the total wave function in "transmitted" and 
"scattered" components. The transmitted wave was evaluated explicitly, but the scattered wave 
was calculated only in the asymptotic limit. The results agreed with those of ||Q. 

Another derivation of the AB scattering amplitude and of the wave function ne ar the forward 
direction has been provided by Dasnieres de Veigy and Ouvry , this time by analyzing the action 
of the propagator on an incide nt plane wave. The procedure resembles those of [|Q and ||, and 
yields identical results. 

One of the few genuine time-dependent analyses of the AB effect has been given by Stelitano 
in a recent work This author decomposes the propagator in positive, negative and zero angular 
momentum components, and calculates each term independently in the asymptotic limit, following 
a procedure similar to the original method of Aharonov and Bohm [jij. His initial wave functio 
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n is a well-localized Gaussian wave packet that approaches the solenoid. The final wave function 
is calculated as the convolution of the propagator with the initial wave packet. The resulting 
scattering amplitude is the same as in the time-independent analyses. The wave packet undergoes 
a self- interference along the forward direction analogous to the one found in GJ and jjj. 

It is clear that there is wide agreement about the results of these calculations. What we present 
in this work is a complete evaluation of the AB wave function based in the results of ||. Our 
method is also applied to calculate the AB propagator. Therefore, both time-dependent and time- 
independent problems can be treated in a similar way. We shall produce analytic results in terms 
of confluent hyper geometric functions which are exact in any angular region, thereby being valid at 
the forward direction. This method was put forward recently in ||], where it was used to calculate 
the wave function in (2+l)-dimensional quantum scattering by a fixed spinless source. It is based 
on Pauli's article on the diffraction of light by a wedge limited by two perfectly reflecting planes 

This article is organized as follows. In Section 2 we shall review Jackiw's solution of the 
Schrodinger equation of a charged particle in presence of a magnetic vortex ||. In a time- 
independent analysis the scattered wave function is given as an integral. The time-dependent 
propagator will be shown to be proportional to the same integral, which is therefore the obje ct to 
calculate. That calculation is carried out in Section 3. The result is applied to the determination 
of the time-independent wave function in Section 4. The time-dependent analysis is presented in 
Section 5. The last Section contains a discussion of the results and our conclusions. 



2 The wave function 

In this Section we shall review briefly the solution given in || to the Schrodinger equation for 
the time-independent AB scattering problem. If we consider that the magnetic vortex coincides 
with the z axis, and that the incoming charged particles approach perpendicularly the vortex, 
the scattering process is esentially two-dimensional. In this situation a possibl e choice of vector 
potential is 

A(r) = — V0 = — uVe, (1) 
2-7T e 

where <3> is the flux carried by the vortex, v stands for the "numerical flux" defined as v = e&/27rhc 
and 9 is the pol ar angle of cylindrical coordinates. The Hamiltonian that defines the dynamics o 
f the system is 

We -- 2 



W =2Ml P " A(r) J (2) 

with A(r) given by Eq. (||). The quantum mechanical evolutio n of the particle- vortex system is 
determined by the time-dependent Schroding er equation, 

ih^(r,0;u;t)=Hy(r,e;v;t). (3) 

In a time-independent scattering process the Schrodinger equation (§) reduces to an energy eigen- 
value problem if time is factorized in the usual way, 

HV E (r,0;v) = EV E (r,9;v). (4) 



2 



Im z 



Rez 



-371/2 



TC/2 



Figure 1: The Schlafli contour 



Since the spectrum of the Hamiltonian TC is continuous, with all energies E > 0, the interest resides 
in the eigenfunctions. These eigenfunctions are found to be || 

\~M~ { [v] 00 1 

V E (r,e ] u) = J— z I £ e - i ^)/ 2 J 1/ _ i (fcr)e^+ £ e in ^ +j ^ 2 Jj^ u (kr) e ije \ . (5) 

V 2 ™ [i=-OQ i = H+i J 

where /c 2 = 2ME/Ti 2 . The brackets [ ] indicate integer part. The sums can be performed with 
the help of the Schlafli representation for the Bessel functions, whose contour of integration C s is 
depicted in Fig. 1 (e is a small positive value for Imz): 

J a (x)=e lan ' 2 [ fkg-fccosjgfaa ^ 

JCs 2vr 

The sums in Eq. (||) are now geometric and can be evaluated in the usual way. After this 
summations and some changes of variable (we refer the reader to || for details), the wave function 
can be represented by 

V 2-kU Jc 27r 1 - e lz 

The new contour C is shown in Fig. 2. It is easy to verify that the wave function ([?]) is single- valued, 
i.e. periodic in 6 wit h 2ir period. The contour avoids the poles z = 0mod(27r). An equivalent 
contour is depicted in Fig. 3. This new contour consists of a loop that encloses the pole z = 0, 
and two vertical lines. This separation corresponds to a decomposition of the total wave function 
^^(r, 9; v) in two components, identified as the transmitted wave (closed contour) and the scattered 
wave (straight lines), 

^ E (r,9;v) 
V tl (r,9;v) 
^ sc (r,9;v) 



M 

— j(*tr(r,0;i/)+ttBc(r,0;i/)) 



V 2irh 

^ikr cos tt) 



oo 



e i[u]e sin(Tw) / j^ fcrcoshi/ ^ . (8) 



7T 
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Figure 2: Contour of integration 
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Figure 3: An equivalent contour 
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The transmitted wave has been evaluated by Cauchy's residue theorem. The scattered wave is 
an integral that can be calculated in the large kr limit by means of a saddle point approximation 
except if 9 = 0. If 9 = the integrand develops a singularity exactly at the saddle point y = 0, 
thus invalidating the procedure. 

Before proceeding to the evaluation of the integral, we demonstrate that it also appears in the 
AB propagator. It has been shown in |8j that this propagator can be written as a Bessel series: 

I e^e-^J^^f), (9) 

v J n=— oo v " 

where <j) = 9 — 9'. The sum can be transformed into an integral by means of the Schlafli represen- 
tation for the Bessel functions, Eq. (|6|). Following the same steps as for the time-independent wave 
function, the AB propagator can be decomposed in two terms that correspond to the transmitted 
and the scattered waves: 

G, c (r,0; r',0';M) = -^exp {i|^(r 2 + e^sin^}^" 1 

1X3 f rarr' 1 e y ^ 

dyexp<i— — coshw> — . (10) 

\ Tit y j l + e»-»* 

— oo 

The computational difficulties lie now in the region <j> ~ ±7r, whe re a straightforward saddle-point 
approximation is ill-defined. The conclusion of this Section is that the evaluation of both the time- 
independent wave function and the time-dependent propagator reduces to a careful calculation of 
the same integral, which we shall denote by I(p, 6, {^}), 

i(pAM) = j^jye^y-^^. (ii) 
3 Calculation of an integral 

In this Section we shall calculate the integral I(p, 9, {^}) defined in Eq. (|Tl|). We are interested in 
an evaluation of this integral that be valid even if the angular variable 9 6 [0, 2n) takes values close 
to or 2tt. 

Let us change to a new variable s = e~ ln ^\/2 sinh |, which extract s a Gaussian factor 
exp(— ps 2 ) in the integrand. The path of integration can be taken as the real axis, so that the 
function I(p, 9, {v}) be 



/oo 
ds 
-oo 



_„«2 1 

e 

oo 



1 — cos 9 + is 2 



cos rj — cos 9 e l ^' n 
x e «v-o) _ i ^T- (12) 



We introduce now the following notation: 

cos rj — cos 9 e 1 ^ 11 
A s >^ = e i(v-e)-i 



a(9) = e _i7r (l-cos0). (13) 
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The choice of phase in the definition of a(ff) is not trivial since thi s function will appear below 
under a square root. The function I(p, 9, {v}) is now 



/(p, 9, {v}) = e- in / 4 V2e ip / ds e~ ps 



f(s,0) 
ia{9) + s 



2 4 



(14) 



The obvious procedure would be to expand the whole integrand in Eq. (|l4|) except the Gaussian 
exponential, in powers of s and evaluate the Gaussian integrals. The result would be ill-defined at 
9 = 0. An alternative procedure is to expand f(s,9) only, according to the following formula: 



f(s,9)=J2^A m (9y 

m=0 



(15) 



This expansion is well-defined for all values of 9. The odd powers of s do not contribute to the 
integration in Eq. (|l4]). The functions A m (9) for m = 0, 2 are 



MO) = 5(^-1) 



A 2 {9) = - (2M - 1) 3 + e t0 - 2{v} + 2{v}e 



J6 



(16) 



Now we insert the expansion fll5| ) in Eq. (|14| ) and consider p 7^ 0. After the change of variable 
s = rp^ 1 / 2 , the function I(p, 9,{u}) reads 



m=0 



Jim 



dre~ T 



(17) 



ia(6»)p + r 2 

This final integration can be performed in terms of confluent hypergeometric functions as follows: 



Fo ( 1, —to H — , ix ) = r ( to 



1 



dr e T 



2 T 



2m 



IX + T 



2 ' 



(18) 



Therefore the final result is 

00 

I(p,9,M) = e^VSe* i m A 2m {0)p^ 



m=0 



xT[m-l)F 2 [l,-m+^,ia(e)p 



(19) 



The behaviour of the confluent hypergeometric functions defined in Eq. ( |18|) for large or small \x 



is 



F 2 (l,-m + §,ix) «(m,-±)(ix) 1 1 - (to + ±) (ix) x + - 
F 2 (l, -m + §,ix) « l + (-m + §) 1 ix + (-m + §) 1 f-m + §) 1 (ix) 2 + 



x » 1 



+vrr (to - i) 1 e -*(W2+ 



|x| « . 
(20) 

By means of these formulas we can evaluate I(p, 9, {^}) in three limiting cases whi ch will be 
important in the next Section: 
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3.1 \pa(9)\ ^oo 
The leading term in this limit is 

I(p, 9, M) = - ,/f e^/V-^ + • • • (21.) 
V 2p sin | 

3.2 p > and 9 = 0+ or 27r~ 

As a shorthand we shall refer to these two possibilities as 9 = + and 9 = 0~, keeping in mind 
that is an angular variable that takes value in the interval [0, 2tt). The leading contributions to 
HP, M) are: 

I(p, ± ,M) = ±™j p + (2M - l)Me^ + • • • (22) 

V 2 p 



The value of the discontinuity I(p, , {z/}) — I(p, , {z/}) could have been determined from (|ll|) 
by means of a simple distributional calculation based on the identity 

V-^—±7ri5(y), (23) 



1 — e~y T it 1 - e-v 

where e — > + and T 7 denotes the principal value. The 0(p -1 / 2 ) term in the right side of Eq. ( p2| 
vanishes if {i^} = 1/2. That this is exact follows from a direct calculation using Eq. (|23]): 



3.3 p = and {u} ^ 



/(p.O^HiW (24) 



This last case can be calculated directly from the definition of I(p,6, {^}) and t he identity (p3|). 
The result is 

1(0, 9, {u}) = -7re i{u}d e- i7IU — ^— . (25) 

sm(7r^J 

To conclude this Section we wish to remark that we have found a unified expression for I(p, 9, {z^}) 
that embraces all particular cases as limiting values o f the parameters 9 and {v}, with the only 
exception of p = 0. We shall show in the next Section how these results can be succesfully applied 
t o the calculation of the AB wave function. 



4 Application to time-independent AB scattering 

We recall the solution to the static AB scattering problem found in [||: 

m tr (r,9;v) = e i^-^) e ikrcose 

% c (r,9;u) = jWsmfrvfr^IfaOM), (26) 

where I(p, 9, {z^}) is the function evaluated in the preceding Section. Therefore we can proceed 
directly to calculate the scattered wave function ^ sc (r,9; v) in som e cases of interest. 
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4.1 kr(l — cos 9) — > oo 

This is the large At limit away from the classical scattering direction 9 = 0. Using Eq. (^) we find 



T / /. % / « f-r i/. .0 + 71"! sin(7T^) 

^sc(r, 0; i/) = W ^ exp + t }> — ^ e + ' ' ' ( 27 ) 



sin 3 



The large-r asymptote of ^ S c( r , 9; u) defines the scattering amplitude f(k, 9) through the formula 
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M, sc (r,9;v) r ^J-f(k,9)e tkr . (28) 



r 



Comparing Eqs. ( P7|) and (28) we extract the usual AB scattering amplitude: 



f{k, 9) = ^= exp U[u] + 1)9 + 4} EH (29) 



IV 2 



sm 2 



Other definitions of the scattering amplitude are possible, depending on what is meant by "incom- 
ing" and "scattered" waves. For instance we may consider that the incoming wave function is a 
plane wave e lkrcose , and define the scattered wave as 



% c (r, 9; v) = VP 8C (r 5 9; v) + * tr (r, 9; v) - e lkrcosy . (30) 

i.e. the total wave function given by Eq. (||), minus the incoming plane wave. As explained in || ; 
the large r behaviour of this new scattered wave defines a scattering amplitude that differs from 
(]29D in a 5(9) -function, thus reproducing Ruijsenaar' s result in Ref. [11]. 



4.2 9 = ± and r > 



We are going to determine the wave function at the forward direction with the help of Eq. (22) 
The scattered wave develops a finite discontinuit y, 



* sc (r, ± ; v) = ±i sin(vrz.) e lkr + {2{v} - 1) W — sin(vri/) e ikr + ■■■ (31) 



As for the transmitted wave, we find directly from Eq. (|26|) that it presents a matching discontinuity: 

^ tr (r,0 + ;^) -* tr (r,27r _ ;z/) = -2i sin(vr^) e ikr . (32) 

Hence the total wave function x I / tr(^ ) 9; v) + ^sc^, 9; v) is continuous at the forward direction as it 
should be. It is possible to show that not only the discontinuities in the scattered and transmitted 
wave functions cancel, but also that all the discontinuities in the derivatives of one of them cancel 
against the discontinuities in the derivatives of the other. The outcome is that the analytic ex- 
pression for the AB wave function presented in Eqs. ( |2"6| ) and ( |T9| ) is smooth even at the forward 
direction, and can be used to examine the behaviour of the wave function for any value of 9. 

It can be said that the decomposition of the total wave function in "scattered" and "transmited" 
components does not make sense at the forward direction (or in general at the classical scatterin 
g directions, see Q for the situation in 2+1 dimensional gravitational scattering). We have seen 
that if we allow for separately discontinuous scattered and transmitted components of the total 
wave function, that decomposition can be sustained and yields a perfectly regular expression for 
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the total wave function. Therefore we can calculate its value at 9 = by averaging between + and 
2-7r~. The result i s 



*(r, 0; v) = cos(Tru)e ikr + (2{v\ - lW — V~ sin(vr^) e ikr + ■■■ (33) 

y lukr 

It is interesting to note that the total time-independent wave function at the forward direction 
vanishes exactly if {v} = 1/2 as a consequence of Eq. (|24|). That value for the fractional part of 
the numerical flux causes maximal scattering outside 9 = 0, and a tot al extinction of the wave 
function at 9 = 0. This can be understood as a consequence of probability conservation. 



4.3 r = and {v} ^ 

Here we shall check that the wave function vanishes at the location of the magnetic vortex. Going 

] we find that the scattered and transmitted wave 



back to Eq. (26) and taking into account Eq. 
functions cancel out: 



*tr(O,0,i/) = e^ 6 -^ 



(34) 



5 Time-dependent scattering 

A similar analysis can be performed for the time-dependent scattering problem . The incoming 
wave function will be now a well-localized Gaussian wave packet of width £ that approaches the 
magnetic vortex from a long distance ro » £ at t = 0, with momentum k : 

\& (r',0';0) = -j^exp ^ikr cos 9' - ^(r' 2 + r% + 2rr' cos 9')^. (35) 

In what follows we shall consider that kr^ » 1 and k » tq/^ 2 . The angle 9' is strongly confined 
at 9' ~ 7r, and similarl y r' ~ ro. These approximations are needed to render feasible the following 
calculations. Therefore we shall not be able to ascertain the behaviour of the outgoing wave packet 
but in an approximate sense. 

Having said that, the outgoing wave packet will be calculated as the convolution of the initial 
wave function (|35|) with the AB propagator, 

*(r, 9; u;t) = J r'dr' J d0' G(r, 9; r', 9'; v\ t)* (r', 0'; 0). (36) 
o o 

We are therefore interested in an explicit evaluation of the propagator found in Sect. 2. The 
"scattered" propagator G sc can be written i n terms of the function I(p, 9, {v}) as 

G sc (r, 9; r' , 9'; u; t) = eX p {^(r 2 + r' 2 ) j sin^}^- 1 / (^, + tt, M^j , (37) 

where <j) = 9 — 9'. Following the same steps as in the previous Section, we can evaluate G sc in the 
following situations: 
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5.1 (f) ~ n or ~ n+ 

This case will be referred to as the "forward direction". What we mean by this needs some 
clarification. If we are considering the propagator alone, there is no doubt that (j) = ±7r does indeed 
correspond to the forward propagation of an incoming particle. However, since our initial wave 
packet has a finite size, the forward direction is not well defined for the outgoing wave packet. The 



difficulty can be traced back to the fact that 9' in Eq. (36) is a variable of integration, and therefore 
the condition <j> = ±7r does not define a fixed value for the outgoing angle 9. A more realistic 
description makes use of the two length scales introduced in the scattering process by the initial 



wav e packet (|35|): £ and rp. Qualitatively, the relevant values of 9' are 

(38) 



arctan + arctan « 



and therefore the condition <f> = it or — ir + means that the outgoing angle 9 is contained in the 
"forward cone" 



arctan — , arctan — 



(39) 



Aside from these considerations, the situation is very similar to the time-independent scattering at 
the forward direction, and similar remarks apply. The propagator G sc splits in two parts, 

G sc (r,e;r',e';u;t) = Gf(r,9;r',9';v;t) + G2(r,9;r\9';v;t) 

Gf(r,6;r',e';v;t) = exp |^(r + r') 2 } simVvr) 

G 2 M;r',0W) = ( 2 M-l)^|^exp{z^(r + r') 2 }sin(^) + ... (40) 

The double sign ± in G\ is + if <f> ~ — ir + and is — if <fi ~ ir~ . We see that there is a finite 
discontinuity in G\ analogous to the one found in the time-independent analysis. This discontinui 
ty is cancelled out by a similar discontinuity in Gt r , which is the propagator associated to the 
transmitted wave: 

G+(r, 9; r', 9'; v- 1) - G t "(r, 9; r', 9'; v- 1) = exp |i^(r + r') 2 | sin(^), (41) 

where the superindices + or — have the same meaning as in Eq. Thus the total propagator is 

continuous at the forward direction. Also, all discontinuities in the derivatives of G sc at the forward 
direction cancel with those of Gtr- The calculation of the total wave function at 9 = can therefore 
be performed by averaging between 9 = + and 9 = 2ir~ , with the following result : 



¥(r, 0; v- 1) = cos(nu) ^ free (r, t) + (2{u} - 1) y — sin(Tr^) V? Ice (r, t) + • • • (42) 

The notation ^f ree (r, t) stands for a free wave packet that propagates along the forward direction, 
not to be confused with ^§ ee (r, t), which propagates radially, 

oo 2ir 

*free(r,t) = J r'dr' J d0' exp j-ifcr' - ^(r' -r ) 2 | 

exp| ? — (r + r) j 



o o 
rn 



2-iriht 
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1 free 



oo 

(r,t) = J dr'exp |— i/cr' — ^2 (/ — ro) 2 | 


exp{i^(r + r') 2 }. (43) 



m 

x. 



A comparison between Eqs. fl33|) and (|4^ ) shows that the behaviour of the time-dependent and 
time-independent wave functions is similar at the forward direction. The leading contribution is 
modulated by th e same factor cos (tti/); also, the ©(r^ 1 / 2 ) -term are identica 1. 

5.2 <p 7^ ±7rmod(27r) and kr — ► oo 

The outgoing angle is now outside the forward cone. In this angular region the transmitted wave 
is subdominant to the scattered wave, so that it suffices to take G sc as propagator. With the help 
of Eq. (|2l| ) we find that the leading contribution to G sc is 



G^r, 0; r-, ff; v, t) _ -^=|=exp + r <) 2 } «**-n«,+r)£^ + - (44) 

An exact evaluation of Eq. (|36j) with G sc and ^ given by Eqs. (44) and fl35|) respectively is not 



known to us. In this circumstance we can resort to the approximation of taking 6' = ir in G sc , but 
not in the wave function v&o- Within this approximation the leading (large kr) contribution to the 
scattered wave function can be easily calculated, with the following result: 



* sc (r, 0;^) = ^^- exp \i ([„] + I) e] *g ee (r, t), (45) 

where \&^ ee (r, t) is the freely propagating one-dimensional (radial) wave packet defined in Eq. (|43|) . 
Again, the AB scattering amplitude can be derived from Eq. (|45|). The result is of course 



Eq. (H). 



6 Conclusions 

The general solution || to the time-independent AB wave function, Eq. (||), has been evaluated as 
a series of confluent hypergeometric functions ([H], p(f ) . The same procedure can be applied to the 
AB propagator, wherefrom the propagation of a Gaussian wave packet has been determined. The 
method is inspired by a recent work || on (2+l)-dimensional gravitational scattering, which in its 
turn was based on Pauli's article [ 10 1 on diffraction of light. We have applied that procedure to 



determine the AB wave function far from and along the forward direction. Our results agree with 
those of previous analyses, bot h in the time-independent or time-dependent approach [||, [|. [|] . 

The factor cos(7r^) present in Eqs. fl33| ) and ( |42| ) is usually interpreted as a self-interference 
between two "halves" of the incoming wave that are split by the magnetic vortex and recombined at 
the forward direction. Each half carries a phase exp(±i-7rz/), whose addition produces the cos(7r^). 
This explanation is only heuristic, sin ce the scattered wave is not taken into account. We have 
seen that the peculiar behaviour of the scattered wave in the forward direction, Eq. ( |3"TD , is crucial 
in rendering the total wave function continuous and single-valued at 9 = 0. 

On the other hand, the separation of the wave function in "scattered" and a "transmitted" 



components becomes arbitrary precisely at the forward direction. The first term in Eq. (31) could 
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as well be considered part of the transmitted wave. Our opinion is that these descriptions of the A 
B effect rely on a semiclassical conception of the wave function, which is supposed to be separable 
in parts that propagate independently. The calculati on presented here recovers the essential unity 
of the wave function after summin g the scattered and transmitted components, and only then we 
find the distincti ve behaviour of the AB wave function in the forward direction. 
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